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Abstract: We investigate a particular version of the Peccei-Quinn (PQ) NMSSM characterized 
by an economical and rigidly hierarchical flavor structure and based on flavored gauge mediation 
and on some considerations inspired by string theory GUTs. In this way we can express the La- 
grangian of the PQ NMSSM through very few parameters. The obtained model is studied numer¬ 
ically and confronted with the most relevant phenomenological constraints. We show that typical 
spectra are for the most part too heavy to be significantly probed at the LHC, but regions of the 
parameter space exist yielding signatures that might possibly be observed during Run II. We also 
calculate the fine tuning of the model. We show that, in spite of the appearance of large scales in 
the superpotential and soft terms, it does not exceed the tuning present in the MSSM for equivalent 
spectra, which is of the order of 10^. 
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1 Introduction 

Low scale supersymmetry (SUSY) is still a good candidate for physics beyond the Standard Model 
(SM) despite a discouraging lack of positive signals in the first run of the LHC. In its uncostrained 
version the Minimal Supersymmetric Standard Model (MSSM) does not provide a solution to the 
origin of the flavor structure of the Yukawa couplings, but hierarchical textures are typical of many 
UV completions, starting from the seminal paper of Lroggatt and Nielsen [1], up to more recent 
developments in L-theory Grand Unified Theories (GUTs), see [2-8] for some early papers. In¬ 
terestingly, L-theory GUTs provide a natural connection to SUSY breaking, as the visible chiral 
matter as well as the messengers responsible for conveying the breaking to the visible sector orig¬ 
inate from the same D7-brane intersection [6]. A consequence of this fact is the requirement that 
matter and messengers share common Yukawa couplings. 

Lrom the phenomenological point of view, a connection between flavor and SUSY breaking 
has been developed in models with flavored gauge mediation [9-12], a branch of gauge mediated 
SUSY breaking (GMSB) [13-19] that links the flavor and messenger sectors of the MSSM. Along 
these lines, one of us recently proposed [20] a very economical model that successfully incorporates 
the advantages of Lroggat-Nielsen-like hierarchical Yukawa structures into the GMSB messenger 
sector. 
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In [20], all the soft terms depend on one extra parameter, h^, which is the coupling between 
the messenger in the 5 representation of SU (5) and chiral matter. Through well known mixing 
effects, the extra coupling can also produce a large third-generation soft trilinear term. At, which 
can be used to enhance the radiative corrections to the Higgs mass. A phenomenological analysis 
of the low-energy limit of this minimal model in the MSSM was presented in [21]. It was shown 
that, despite quite restrictive constraints on the nature of the matter-messenger couplings, which 
arise from the flavor structure, it is easy to find regions of the parameter space characterized by 
particle spectra consistent with the bounds from direct searches for squarks and gluinos at the 
LHC, and with the measured value of the Higgs boson mass, ~ 125 GeV, as well as a number 
of constraints from flavor-changing neutral current (FCNC) in low-energy processes. 

However, the matter/messenger system of Ref. [20] is fully consistent with hierarchical flavor 
structures when the ratio of the Higgs doublets’ vacuum expectation values (vev’s), tan f3, is of 
the order of a few and not larger. As is well known, in the MSSM small values of tan f3 reduce 
the size of the tree-level Higgs mass, which therefore requires substantial radiative corrections 
that imply large soft masses and a consequently high level of fine tuning. Thus, in this paper we 
analyze whether it is possible to obtain a low-energy limit of the model presented in [20] in the 
Next-to-Minimal Supersymmetric Standard Model (NMSSM). As is well known, this is achieved 
by introducing one additional gauge singlet chiral field, S, that couples to the Higgs sector in the 
superpotential. Because the tree-level value of the Higgs mass for low tan /3 values can be easily 
made larger than in the MSSM, the NMSSM seems to be a more natural choice as a low-energy 
effective theory originating from the model in [20] . 

In particular, we investigate here the Peccei-Quinn (PQ) version of the NMSSM, which is 
characterized by the vanishing of the coupling, k = 0 [22-28]. We will build the PQ NMSSM 
effective action based on arguments from its possible UV completion. It is interesting that similar 
models appeared in F-theory constructions [6], so that we will often invoke F-theoretic arguments 
to justify our assumptions. 

The purpose of the paper is twofold. We construct a version of the PQ NMSSM, based on a 
string inspired UV completion and on flavored GMSB, which yields a very predictive framework. 
The effective action is defined by just a handful of free parameters and we show that the obtained 
spectra are consistent with the basic phenomenological constraints, although we also show that 
they are for the most part too heavy to be significantly explored at the LHC. 

We further show that the PQ NMSSM analyzed here requires the unavoidable introduction of 
nonstandard tadpole superpotential and Lagrangian terms, which are linear in the S field, and terms 
quadratic in S as well, all of them characterized by typical scales that are by orders of magnitude 
higher than the SUSY breaking scale. Unlike the standard GMSB soft terms, these new terms break 
the global symmetry of the PQ NMSSM, C/(1)pq, whose existence is at the origin of the k = 0 
choice. In the string theory setup, C/(1)pq is a remnant of a local gauge symmetry which has been 
spontaneously broken. As is well known, the corresponding gauge boson acquires a large (GUT 
scale) mass through the Green-Schwarz mechanism [29]. Here we break [/(1)pq through the vev 
of a scalar field X, which will also play the additional role of the SUSY breaking spurion. 

On the other hand, we show that despite the presence of these extra large-scale terms the model 
does not present fine-tuning levels higher than those already present in its MSSM version, i.e., 
~ 10“^. This is due to a cancellation among some specific terms entering the fine-tuning measure. 
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because of relations determined by our UV completion. Thus, this model is partially immune from 
the mild “tadpole problem” [30-37] of the General NMSSM [38], according to which the extra 
terms dramatically amplify the effective theory’s sensitivity to the high energy physics. 

The paper is organized as follows. In Sec. 2 we introduce the model. We start from the 
flavored GMSB structure of the superpotential and progress subsequently to introducing one by 
one the additional terms that will define our version of the PQ NMSSM. In Sec. 3 we perform a 
numerical analysis of the model. We identify the regions of the parameter space consistent with the 
constraints from the Higgs mass measurement and LHC searches, and we provide some benchmark 
points useful for the discussion and also possible collider signatures. In Sec. 4 we discuss the fine 
tuning of the model, and prove that it is not larger than the present level found in the MSSM, despite 
the presence of tadpole and quadratic terms in the singlet field. We finally provide our summary 
and concluding remarks in Sec. 5. 

Additionally, we provide four appendices to the text, dedicated, respectively, to the GMSB 
calculation of the soft masses in the model; to the explicit calculation of the U{1)pq breaking 
effective terms of the superpotential and soft Lagrangian; to the explicit estimate of the size of our 
constants; and to the explicit form of the tree-level Higgs mass in the PQ NMSSM. 

2 The model 

2.1 Flavored GMSB with hierarchical Yukawa couplings 

A simple model combining the advantages of F-theory model building and flavored GMSB was 
introduced in [20]. Since the visible matter and the messengers have the same origin they should 
also present a common hierarchical structure of their couplings, e.g., a structure of the Froggatt- 
Nielsen type [1]. This reasoning basically implies that the messengers are in 5 representations of 
SU (5). The simplest construction of this type contains three chiral families of 10s, four 5s and one 
extra multiplet 5y, which is necessary to form a vector pair of messengers and cancel the resulting 
anomalies. 

The relevant superpotential terms using SU (5) representations at the GUT scale are 

W D ^ 10,10,(5/r)2 + J^yijlO*5j(5H)2 + (^aj5j)5yX, (2.1) 

i^ij J 

where in Eq. (2.1) flavor indices run as i,j = 1,2,3 and J = 1,2, 3,4, and the subscript “2” 
attached to the Higgs fields highlights that we take into account only the doublet part of the 5h and 
5h Higgs multiplets. In this regard, we recall here that in F-theory GUTs doublet-triplet splitting 
can be obtained through an appropriate choice of internal fluxes (see, e.g.. Secs. 10.3 and 12 of 
Ref. [5]). 

We assume that all the couplings in Eq. (2.1) have a hierarchical structure, i.e., ?/“■« ylt+ij > 
yfj ^ yfj+i F j (and similarly for the yij), and also aj <C oj+i. One can claim [7] that the 

couplings with top-most indices, 2 / 33 , ^ 34 , and 04 are all of the same order. The measured value 
of the top quark mass sets yt = 2/33 to be of 0(1) if the renormalization group (RG) flow does not 
change yt drastically. Note also that, as was mentioned in Sec. 1, the hierarchy of the couplings is 
better enforced for values of tan fi not much larger than 1 , as ^ 33/^34 ~ tan f3 mij/mt. 
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The spurion X acquires a vacuum expectation value (vev), {X), which gives a mass, My = 
{X), to the messengers 5y and 5y, where 5y = The other 5s remain massless. The 

spurion also triggers SUSY breaking through its nontrivial F-term: X —)• {X) + 6‘^Fx ■ We 
assume here that the dynamics of the spurion is governed by the physics of the high scale, so that 
it is effectively a non-dynamical field below the messenger scale. 

After diagonalizing the mass matrices in Eq. (2.1), we can write down the terms with 0(1) 
couplings that will be relevant for our analysis: 

W D ^* 103103 ( 5^)2 + /l3l035y(5/r)2 + 5y5yX , (2.2) 

where y* « 1, (13 = ^34 ~ 1. By integrating out the messengers one can obtain the soft terms of 
the MSSM explicitly displayed in Appendix A. 

2.2 NMSSM GMSB with PQ symmetry 

In the NMSSM one introduces an extra chiral superfield S, whose vev generafes fhe // ferm af fhe 
elecfroweak (EW) scale. We consider here fhe hypofhesis fhaf fhe fields of fhe model are charged 
under a PQ symmefry, C/(1)pq, which infroduces addifional resfricfions on fhe possible couplings. 
In parficular, as we shall see below, if forbids any superpofenfial ferm for S alone (k = 0). We 
also assume fhaf fhe PQ charges of fhe X and S fields are opposite in sign, q(S) = —q(X) (as 
considered, e.g., in [6]), which leads fo cerfain necessary ferms in fhe Eagrangian. Incidenfally, 
t/(l)pQ also helps prevenf fasl baryon decay due fo GUT processes [39]. 

The resulting C/(1)pq charges are summarized in fhe following fable: 



10l,2,3 

5i,2,3, 5y 

5/r 

5h 

5y 

5 

X 

Q 

+1/4 

+1/4 

-1/2 

-1/2 

+3/4 

+1 

-1 


where fhe charges are normalized in such a way fhaf fhe largesl one equals 1. Note fhaf fhis implies 
fhaf fhe spurion X does nof couple fo fhe Higgs fields, while S cannof couple fo fhe messengers in 
dimension 4 operators. Moreover, as menfioned above, fhere cannof exisf polynomial couplings in 
S and X alone so fhaf, e.g., k = 0 in fhe nofafion of fhe usual Z^-symmefric NMSSM [40]. 

In summary, fhe only renormalizable superpofenfial ferm preserving C/(1)pq and confaining S 
is^ 

IU 5 = XSHuHd , (2.4) 

which after SUSY breaking gives rise fo a frilinear soff Eagrangian confribufion, A\, and a scalar 
sofl mass for fhe singlef, m|, 

- L5,soft = XAxHuHdS + ml\Sf , (2.5) 

whose explicif expression in terms of the GMSB parameters is given in Appendix A. 

As one can easily see, U (1)pq is anomalous. It is known that the anomaly can be canceled by 
the Green-Schwartz mechanism [29], which gives a mass to the U{1)pq gauge boson of the order 
of Mqut ■ 

Note that (X) breaks U (1)pq, so that at low energies there appear additional effective C/(1)pq 
breaking interaction terms. We present them in the next subsection and discuss them in detail in 
Appendix B. 

*We drop here the mass term XS, which is unnatural in F-theoretic constructions, see Appendix B. 
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2.3 Effective interactions below the messenger scale 


In addition to (2.4) and (2.5), other contributions to the superpotential and soft Lagrangian might 
arise below My- As we explain in Appendix B, these can be due, for example, to instanton effects, 
exchange of heavy chiral U (1)pq neutral superfields, exchange of heavy gauge bosons, or Giudice- 
Masiero [41] effective terms in the Kahler potential. The most important terms for the low-energy 
dynamics will be those containing powers of the spurion superfield, X, because X receives a large 
vev. 

Additional terms to the superpotential and soft SUSY-breaking Lagrangian take the form 

6Ws = ^fS + , (2.6) 

soft = ~^^sS + h. c.^ , (2.7) 

where we have used the standard notation of [40] and we do not explicitly distinguish between the 
superfields and their scalar components. Importantly, to a very good approximation the relations 

K « 0 and Ak ~ 0 (2.8) 

still hold, even while there are other massive couplings that effectively break (7(1)pq. 

The superpotential and Lagrangian elements in Eqs. (2.6) and (2.7) are best expressed in terms 
of X = My /Mgut and A = Fx/My. They are given by (see Appendix B for full details) 

Sml = -|A^p|Apx^, 

^ j ! = — A|,MguT 

= -2/i,3AMGUTa^^ = -2 ^fA, 

m'l = 2 A|AMgut = -2/i'A, (2.9) 

where /i ^3 is a dimensionless instanton contribution, and A>p is a 0(1) “Yukawa” coupling between 
the X and S superfields and a heavy chiral (/(1)pq neutral superfield, 'k, that is integrated out below 
Mgut ■ Without loss of generality, in what follows we fix Mgut = 2 x IQi® GeV and A^ = 1. 

All the soft terms and superpotential parameters in Eqs. (2.9) are defined at My and then 
renormalized to lower energy through RG flow. The model is therefore described entirely by 5 free 
parameters defined at the messenger scale: 

My (or x), A, h3,X,^F , (2.10) 


where we have used the more familiar superpotential tadpole in place of the numerical instanton 
coefficient Ii 3 . As will be clear below, from a phenomenological point of view one is interested 
in a somewhat limited range of A, so that the number of really free parameters can then be reduced 
to 4, of which only 2, x and set the overall scale of the additional terms in Eqs. (2.6)-(2.7). 

We want to highlight here that 


_ m'^ 

ip 


( 2 . 11 ) 
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which holds up to leading terms (see Appendix C). Equation (2.11) is RG invariant after sup¬ 
pressing suhdominant terms in the RG equations. This characteristic will turn out to he essential 
to reducing the fine tuning of the model hy approximately two orders of magnitude, thus largely 
ameliorating the tadpole problem, as we discuss in detail in Sec. 4. 

3 Phenomenological analysis 

In this sections we present the results of our numerical analysis, which determines the phenomeno¬ 
logical constraints on the parameter space of the model and possible signatures at present and future 
experiments. In practice, the most important constraints come from the measurement of the Higgs 
boson mass, rrih^ ^ 125 GeV [42, 43] and the bounds from direct searches for SUSY at the LHC, 
of which the most relevant are the ones on the gluino mass, rrig > 1.5 x 10^ GeV [44-46], and on 
heavy stable charged particles [47]. 

We calculate spectra and NMSSM parameters with NMSSMTools v4.5.1 [48], which we 
modified to incorporate the / 13 -dependent and A-dependent contributions to the soft terms given in 
Eqs. (A.l) of Appendix A, and the additional masses and tadpoles defined in Eqs. (2.9). 

To guide the scanning procedure we use MultiNest [49]. We scan with flat priors in the 
parameters A, My, X{Ms\jsy), ^ 3 . and tan (3. The fundamental parameter \{My) has been traded 
for the renormalized value of A at the scale of the geometrical average of the stop masses, Msusy; 
the fundamental parameter has been traded for the ratio of the Higgs doublets’ vevs, tan (3, 
through the electroweak symmetry breaking (EWSB) conditions. 

We scan in the following ranges: 

10® GeV < A < 10® GeV, 10^ GeV < My < 10^^ GeV, 

0.1</i3<l, 0.1 < A(Msusy) < 2, 1< tan/3 <8. (3.1) 

As was mentioned above, the LHC lower bound on the gluino mass affects the choice of 
scanning range for A, which we constrain to a narrow interval, as the soft masses of the gluino and 
the bino strongly depend on this parameter. Below A ss 10® GeV the gluino tends to become too 
light, while for A > 10® GeV the model becomes substantially and uncomfortably more fine-tuned. 
The lower bound on My comes from requiring the two-loop GMSB expression of Eqs. (A.l) to be 
the dominant contribution to the soft masses, while the upper bound comes from the requirement of 
being well within the range of validity of the effective theory below Mqut : for My > 10^“^ GeV, 
the first-order expansion at the origin of Eqs. (2.9) is not well-justified. Note finally, that we allow 
for a wide range of A values, to counterbalance our choice of low tan /3. 

In addition to the fundamental parameters we scan in some of the SM nuisance parameters: 
the strong coupling constant, the MS value of the bottom quark mass, and the top quark pole mass, 
which we include in the likelihood function. Eor these we adopt normal distributions based on 
the most recent PDG [50] central values and experimental uncertainties. The scans are driven by 
a Gaussian likelihood function for the experimental measurement of the Higgs mass, where we 
added to the experimental uncertainties, in quadrature, a theoretical uncertainty of approximately 
3 GeV. This choice is motivated by the large uncertainty still present in the calculation of the Higgs 
mass, which includes the choice of renormalization scheme, missing higher order contributions, 
and numerical differences between various existing public codes (see [51] for a recent discussion). 
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Figure 1: (a) The lightest Higgs mass value for the points within the la theoretical uncertainty, 
122 GeV < nih^ < 128 GeV, in the (A(Msusy), tan P) plane, (b) A zoom of the region of large 
A, where we plot the tree-level value of the Higgs mass in the (A(Msusy), tan/3) plane. 


In Fig. 1(a) we show the value of the lightest Higgs mass for the points within the la theo¬ 
retical uncertainty, 122 GeV < rrih^ < 128 GeV, in the (A(Msusy)j tan (3) plane. One can easily 
recognize two different regions. The first is characterized by A Ri 0.1 — 0.5 and tan/3 > 5, and 
the scan seems to give there a Higgs mass slightly low, mh^ < 124 GeV, albeit well within the 
adopted theoretical uncertainty. This “small-A” region shows solutions similar to the MSSM limit 
of the model [21]. The Higgs mass decreases for lower values of tan/3, as a consequence of the 
rapid drop in the tree-level value, but it also decreases slowly for A > 0.4 — 0.5 when tan/3 > 7. 
For the majority of the points this is due to the trilinear term \At \ becoming smaller as A increases,^ 
as the scan cannot compensate this effect by raising Msusy because of our narrow range in A. 

A second region, more interesting for the purposes of this paper, can be found instead at 
A > 0.7 and tan/3 < 3. As Fig. 1(a) shows, the Higgs mass there can comfortably reach the 
measured value and above. However, not all the points in this region are equivalent to each other, 
as an analysis of the tree-level Higgs mass can show. 

In Fig. 1(b) we present a zoomed-in detail of the “large-A” region, where we plot in the third 
dimension the tree-level value of the Higgs mass, whose explicit expression is presented in Ap¬ 
pendix D. Figure 1(b) shows two distinct sets of points, characterized by different properties. For 
A > 0.8 there are points characterized by tree-level masses in the range ~ 70 — 90 GeV. These 
solutions are somewhat similar to models of ASUSY [52], characterized by values of A that can 
become non perturbative before reaching the GUT scale.^ The tree-level mass is enhanced with 
respect to the MSSM value, and as a consequence the needed radiative corrections are less substan¬ 
tial than in the MSSM for equivalent tan /3. Note, however, that we could not find a single region 

^Recall the RGBs for At'. dAt/dt ~ \^Ax [40], where A\ < 0, see, e.g., the second to last of (A.l) in Appendix A. 

^Since there are no direct couplings between the scalar S and the messengers, and S is also a gauge singlet, at one 
loop no additional terms arise in the fix above My ■ The gauge couplings g\ and g 2 enter fix and in the presence of the 
messengers can get renormalized slightly more strongly. We have checked numerically that the effect is however very 
small and does not change the fact that A > 0.8 generally becomes nonperturbative at the GUT scale. 
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Benchmark 

BPl 

BP2 

BP3 

BP4 

Model parameters (at My) 





A 

2.98 X 10® GeV 

5.79 X 10® GeV 

5.92 X 10® GeV 

9.15 X 10® GeV 

My 

2.49 X 10^^ GeV 

5.99 X 10®® GeV 

9.56 X 10®® GeV 

2.48 X 10®® GeV 

hs 

0.92 

0.65 

0.73 

0.23 

X 

1.72 

1.69 

0.14 

0.35 


2.14 X lO^^GeV^ 

2.03 X 10®® GeV® 

1.23 X 10®® GeV® 

3.70 X 10®" GeV® 

Relevant for EWSB (at Msusy) 





tan P 

1.70 

1.40 

7.52 

7.46 


4.33 X 10^ GeV^ 

4.39 X 10'® GeV® 

8.95 X 10® GeV® 

-8.13 X 10° GeV® 


-2.24 X lO^'GeV^ 

-3.19 X 10'® GeV® 

-1.21 X 10® GeV® 

-1.51 X 10® GeV® 

A 

0.82 

0.70 

0.13 

0.33 

Meff 

7639 GeV 

10746 GeV 

3596 GeV 

4048 GeV 


-1746 GeV 

-1353 GeV 

-7635 GeV 

-1075GeV 

m| 

-4.90 X lO’’ GeV^ 

-6.94 X 10'® GeV® 

-1.71 X 10° GeV® 

-2.28 X 10° GeV® 

F 

-1.69 X 10® GeV 

-9.66 X 10®° GeV 

-4.53 X 10® GeV 

-2.89 X 10®° GeV 

m'i 

1.00 X 10^^ GeV^ 

1.11 X 10®'® GeV® 

5.37 X 10®® GeV® 

5.28 X 10®° GeV® 

?F 

1.58 X 10^“ GeV^ 

1.49 X 10®® GeV® 

1.23 X 10®® GeV® 

3.58 X 10®" GeV® 


-9.17 X 10^® GeV® 

-1.70 X 10®® GeV® 

-1.44 X 10®® GeV® 

-6.55 X 10®° GeV® 

Max fine tuning 

~ 10^(A,^f) 

~ 10° (A) 

~ 10" (Cf) 

~ 10" (Cf) 

Spectrum 





,tiee 

74.1 GeV 

27.7 GeV 

84.6 GeV 

84.1 GeV 

ruh^ 

123.5 GeV 

123.0 GeV 

122.6 GeV 

122.3 GeV 

1 

5.86 X 10"* GeV 

2.23 X 10® GeV 

7.11 X 10" GeV 

5.91 X 10° GeV 

nig 

2170GeV 

3959 GeV 

4006 GeV 

5970 GeV 


2311GeV 

3759 GeV 

3436 GeV 

5509 GeV 

m^o 

422 GeV 

813 GeV 

819 GeV 

1266GeV 

m ± 

Xi 

832 GeV 

1563 GeV 

1550 GeV 

2355 GeV 


920 GeV 

2170GeV 

2391GeV 

261 GeV 


Table 1: The model fundamental parameters at My, the corresponding parameter values at Msusy> 
fine tuning, and spectra of four benchmark points discussed in Sec. 3. The parameters in parenthe¬ 
ses give the maximal contribution to the fine-tuning measure. 


of the parameter space in which > M|, so that significant radiative corrections are always 

necessary in this model to obtain the correct Higgs mass. 

Still, these points present somewhat lighter spectra than in the rest of the parameter space, 
although for the most part still too heavy to be significantly probed in Run II at the LHC [53]. The 
lightest gluino mass found by the scan is around 2 TeV, as can be seen in Table 1, where we present 
the parameters and spectral properties of a typical point (BPl). 

Several points in Fig. 1(b) are characterized by a relatively low value of the tree-level Higgs 
mass, ss 20 — 60 GeV, despite sizeable values of A. This can be understood from Eq. (D.3) 

in Appendix D, which shows that for low tan/3 the MSSM-like part of the tree-level Higgs mass 
is augmented in the PQ NMSSM by an extra term which is approximately Sm\_^ e, 

where e = assumes in our model values in the range ~ 0 — 0.8. Thus, the tadpole 

ratio significantly affects the value of the tree-level mass: if ^5 <C CfI^' the tree-level mass is 


































approximately equal to its typical MSSM value, as is exactly the case for these points. 

For A > 0.7, these solutions can still quite comfortably reach rrih^ ss 125 GeV thanks to very 
large A-dependent radiative corrections. A representative point of this region (BP2) can be found 
in Table 1 . One can see that BP2 is characterized by a large messenger scale. My ps 10^^ GeV, 
which is necessary to lift the heavy Higgs masses up and boost the Higgs-loops corrections to the 
lightest Higgs mass. As can be expected, these points are marred by extremely high levels of fine 
tuning, because of the large sensitivity to A-driven radiative corrections. Thus, we will not discuss 
them further in what follows. 

The third benchmark point in Table 1 (BP3) is representative of the small-A, almost MSSM- 
like region shown in Fig. 1(a) and discussed above. With respect to the large-A region (see BPl), 
the lightest points in the small-A region tend to present gluino and neutralino masses approximately 
a factor 2 heavier than the lightest points with large A. Note that the fine funing of fhe small- and 
large-A regions are very comparable. 

A discussion of fhe fine funing of fhe model is presented in Sec. 4. The model is funed af 
teas! fo fhe level of 5 x 10^ - 10^, buf if is interesting fo nofe fhaf, parficularly for fhe poinfs of 
fhe small-A region, which presenf equivalenf MSSM solutions in fhe limif A —)> 0, fhe fine funing 
is nof higher fhan if would be in fhe MSSM, even if in fhe PQ NMSSM presenfed here we had 
fo infroduce addifional large scale terms fhaf enter fhe EWSB condifions: ^ 5 , /r', and m'g. 

In facf, we will show in Sec. 4 fhaf UV relations befween fhese paramefers allow for significanf 
cancellations in fhe fine-funing measure, so fhaf fhis does nof increase proportionally fo fhe fypical 
scale of fhe exfra terms. 

The paramefer ( 13 , which is responsible for fhe mixing befween fhe messenger and maffer 
sector, cannof span fhe entire range of Eq. (3.1): fhere is a gap in ifs allowed values, which was 
already observed for fhe MSSM in [21]. The reason can be easily understood by a quick analysis 
of Eqs. (A. 1) in Appendix A. The sfaus can become fachyonic when a cancellation lakes place 
befween terms dominated by the gauge couplings and those dominated by / 13 . The position and 
range of this gap strongly depend on the values of A, and the gap tends to become smaller for larger 
A values. 

In Eig. 2(a) we show the distribution of the stau mass as a function of for the small-A 
(magenta diamonds) and large-A (indigo triangles) regions. The stau mass can drop drastically as 
one approaches the critical values of giving rise to this cancellation, to the point that it becomes 
the Next-to-lightest SUSY particle (the ESP is the gravitino), as a comparison with the neutralino 
mass, plotted here with gray circles, shows. 

This brings about the interesting possibility of testing selected /13 regions by searches for heavy 
stable charged particles. We plot in Pig. 2(a) with a black dashed line the 95% C.E. lower bound 
on the mass of the stable stau from an analysis of long time-of-flight (TOP) to the outer muon 
system and anomalously high (or low) energy deposition in the inner tracker at CMS, with 8 TeV 
data [47]. The search excludes cross sections of the order of 1.5 fb or above for pair production of 
~ 300 GeV stable staus, which is in agreement with what one obtains in this model with typical 
spectra and a stau NESP with a mass of 300 GeV. We thus expect that the corresponding search in 
Run II will start biting into specific regions of our paramefer space: /13 « 0.2 or 0.5 for A < 0.4; 
(13 0.4 for A Ri 0.85 — 0.9. We show a represenfafive poinf wifh sfau NESP (BP4) in Table 1. 

We conclude fhis secfion by comparing in Pig. 2(b) fhe projecfed [53] exclusion reach for 
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(a) (b) 

Figure 2: (a) The distributions of the lightest stau mass, mf^ as a function of for the points of 
the small-A region (magenta diamonds) and of the large-A region (indigo triangles). The neutralino 
mass is plotted in the background, as gray circles. The dashed black line shows the 95% C.L. lower 
mass bound on heavy stable charged particles from the tracker + TOF analysis at CMS 8 TeV [47]. 
(b) The ATLAS projections for gluino searches at the LHC 14 TeV [53] compared to the viable 
points of our model. The dot-dashed black line shown the 95% C.L. expected exclusion reach with 
300 fb“^ and the dotted red line the reach with 3000 fb“^. 

gluino masses at ATLAS 14 TeV with our model’s predictions. Again, points of the large-A region 
are shown as indigo triangles and points of small A as magenta diamonds. The 95% C.L. expected 
reach with 300 fb“^ of integrated luminosity in searches with jets -i- missing Et is shown as a dot- 
dashed black line. As is well known, in GMSB models specific LHC phenomenology is strongly 
affected by the effective coupling of the neutralino (which is predominantly bino-like in our model) 
and the gravitino. However, gluino searches at 8 TeV have shown that the mass bounds do not 
significantly depend on whether the neutralino is long-lived enough to escape the detector [44, 46], 
or it decays promptly producing photons [45]. Thus, for the purpose of this paper we accept the 
projected reach of [53] as a reasonable estimate. 

One can see in Fig. 2(b), that the points of the large-A region might begin to be probed at 
the very end of the LHC Run II, particularly if the High-Luminosity LHC [54], whose reach [53] 
is shown as a dotted red line, is approved. However, Fig. 2(b) more realistically shows that the 
structure devised here yields spectra too heavy to be significantly tested at the LHC, and a high- 
energy collider will be necessary. It might be interesting to study the reach of a 100 TeV collider 
for this model, and we leave this for future work. 

4 Fine tuning 

We dedicate this section to the calculation of the fine tuning of the model. As was mentioned 
in Sec. 3, it turns out that this is smaller than what one could expect from naive estimates based 
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on a rough order-of-magnitude analysis. We take the maximum of the fine tuning due to each of 
the fundamental parameters of Eq. (2.10), which we calculate according to the Barbieri-Giudice 
measure [55, 56]. 

By indicating the parameters of Eq. (2.10) collectively as pi, one must calculate the max of the 


d log M| 
<9 log Pi 


Pi dMy 

— 7) (order of magnitude only). 

Mz ^Pi 


(4.1) 


Given the three EWSB conditions that determine the values of the Higgs vevs,"^ 


EWSBj(M|, tan^,s,/r',^ 5 , m'|...;pi) = 0 forj = 1,2,3, (4.2) 


one can apply the chain rule, 

dPk _ r9(EWSBj) 
dpi dPk 

where the 3x3 matrix in the jk indices is built out of the partial derivatives with respect to 
Pk = {Mz, tan (3, s} , and one is interested in the numerical value of the first row in the system 
(4.3), dMljdpi. 

Despite the presence of large contributions to EWSB from the tadpoles and other terms, a 
numerical calculation shows that the partial derivatives are sensitive to cancellations. To see this, 
we parametrize two large parameters, ^5 and m'g , in the following way: 

Cs = -aA^F, m's = -hAp! , (4.4) 


a(EWSBj 

dpi 


(4.3) 


where a and b are arbitrary real numbers, and we proceed to the calculation of and Amy for 
changing a and b. My, which directly enters the definition of //', and are the largest remaining 
free parameters in the theory. 

We show in Eig. 3 a plot of logj^Q in the (a, b) plane for BPl. A plot of log^g Amy shows 
similar behavior and we do not present it here. One can clearly see that the fine funing is reduced 
by more fhan fwo orders of magnifude when a Fib. 

One can see from Eqs. (B. 8 ) and (B.9) of Appendix B fhaf in our model 6 = 2 is a consequence 
of breaking SUSY via fhe spurion. As a consequence, according fo Eig. 3, one needs o ss 2 fo lower 
fhe fine funing fo ~ 10^. As is explained in Appendix B after Eq. (B.4), we make here a reasonable 
assumpfion, fhaf fhere exisf selection rules in fhe sfringy UV completion fhaf allow one fo generate 
only fhe insfanfon confribufions fhaf lead fo a = 2. Note fhaf, as one can see in Eqs. (B.13) and 
following paragraph, fhe ferms depending on are much smaller fhan fhe insfanfon confribufion 
fo ^ 5 , so fhaf 7 ^ 0 will nol cause significanl deviafions from o = 2 . 

The cancellafions described here have fheir origin in fhe specific form of fhe free-level EWSB 
conditions. They lead fo fhe fine-luning levels of fhe benchmark poinfs in Table 1, and in general 
similar values are obfained over fhe whole viable paramefer space. Thus, as was anticipated in 
Sec. 1, our model is parfially protected from fhe fadpole problem, provided an appropriate relation 
of fhe kind of (4.4) is given by fhe UV completion. 

'^The explicit forms of EWSBi, EWSB 2 , and EWSB 3 can be found in Eqs. (C.1)-(C.3) of Appendix C. 
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Figure 3: A plot of log^Q for BPl of Table 1 as a function of the arbitrary real coefficient o, b. 

5 Summary and conclusions 

In this paper, motivated by the desire of maintaining a Yukawa flavor structure in agreement with 
UV completions based on F-theory, we have used flavored gauge mediation to construct a specific 
version of the Peccei-Quinn NMSSM. Considerations concerning the stringy UV completion lead 
to a very predictive version of the model, which depends on a few free parameters. 

We have performed a thorough numerical study of the parameter space of the model and 
confronted it with the phenomenological constraints. Our findings are supported by analytical 
arguments. We showed that the PQ NMSSM exhibits unusual properties, specifically it requires 
relatively large nonstandard terms: tadpole coefficients ^ 5 , and quadratic terms depending on 
/x' and m'g. In spite of this, we showed that the fine tuning is not greater than in the MSSM thanks 
to a special relation between the ratios ^s/^F and m'g/n', which originates in the specific UV 
completion considered here, and which results in a cancellation in the fine-tuning measure. Thus, 
we showed that this version of the the PQ NMSSM is partially immune from the mild tadpole 
problem typical of many general versions of the NMSSM. 

We found that the model can easily accommodate the Higgs mass at 125 GeV, but the spectra 
can possibly be in reach of the LHC only in a region of the parameter space characterized by 
A > 0.8 which, as is well known, implies that the renormalized value of A at the GUT scale might 
incur a Landau pole. In general, we were not able to find a single parameter space region over 
which the tree-level value of the lightest Higgs mass is larger than the Z boson mass, although it 
can be enhanced with respect to its MSSM counterpart in regions of large A and small tan /3. In 
all cases, however, ruh^ « 125 GeV requires substantial radiative correction and a certain level of 
fine tuning, of the order of 10 “^, cannot be avoided. 

In spite of these interesting formal properties, we found that the phenomenology of the model 
is quite standard, with typical spectra being too heavy to be significantly probed at the LHC. How¬ 
ever, small regions of the parameter space exist featuring gluino masses around 2 TeV, for a bino- 
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like neutralino mass of ~ 400 GeV. These might begin to be probed at the end of Run II in gluino 
searches. For particular choices of the parameter / 13 , which in our model determines the mixing 
between the messengers and matter sector, the NLSP is a stau with a mass > 300 GeV, and a cross 
section right in the ballpark of the range currently probed by searches for heavy stable charged 
particles at CMS, which will be able to further constrain part of the parameter space with 13 and 
14 TeV data. 

This work can be extended in several directions. An open question remains concerning the 
dynamics of the spurion X, which should result from a UV completion of the model (supergravity 
or a string compactification with fluxes). It is of course possible that most of the necessary ingre¬ 
dients are contained in Appendix B, with the exception of some self-interaction terms in X. In that 
case, including spurion quanta into the analysis should not pose a difficult task. 

Another possibility would be to enlarge the size of the corrections given by exchange of the 
C/(l)pQ gauge bosons by increasing the strength of the gpg coupling. This would amend the 
soft terms given in Appendix A by possibly substantial corrections. Finally, one could lower the 
messenger mass, which would result in a modification of the soft masses by one-loop terms. 

The list of modifications is of course much longer, and we leave this extended discussion to a 
future publication. 

Acknowledgments 

We would like to thank Tomasz Jelinski and Andrew Williams for helpful discussions and inputs. 
K.K. is supported in part by the EU and MSHE Grant No. POIG.02.03.00-00-013/09. K.K. and 
E.M.S. are funded in part by the Welcome Programme of the Eoundation for Polish Science. The 
use of the CIS computer cluster at the National Centre for Nuclear Research is gratefully acknowl¬ 
edged. 


- 13 - 


A Soft terms 


At the leading order, all soft masses and trilinear couplings are generated at the messenger scale 
My = {X). While deriving the explicit formulas, we will make some simplifying assumptions. 

First of all recall that Eq. (2.2) is defined at Mgut- Evolution to My leads to further mixing 
between matter and messengers. However, we have explicitly checked that this has a very mild 
influence on the structure and values of the SM Yukawa couplings and /13 at My . Secondly, we 
take into account only the leading superpotential contributions, Eq. (2.2), thus disregarding all other 
contributions from Eq. (2.1), which are usually much smaller. In particular, we neglect the terms 
proportional to the bottom and tau Yukawa couplings. 

In the end, all the soft terms depend only on A = Fx/{X), A, /13 = ^ 34 , yt = 2 / 33 , and 
the gauge coupling constants, while the dependence on the scale My appears through the RG 
evolution. 

We adopt here the general formulas derived in [11, 57] in the presence of MSSM-messenger 
superpotential interactions, which provide a good approximation for the case at hand. The resulting 
soft terms are thus 




'Q>i3 
2 


A 2 

{a\ + Ahal + SOag) + [l05|/r3p - dvr (7ai + 45a2 + SOas)] } 

^SijTr'^ {oy + 5al) - , 


{'>^3)1 = 


SOtt^ 

A 2 


^[5,vr^K + 20a^)], 

{aj + 5ai)] , 

A 2 

= f5407r4 [35|/i3p - 127r (Sai + 5a2)] } , 

3A2 








leovr^ 

3A2 


TT^ial + 5al) - ^|/i3p|ytP - 

0 D 


A, = 


Ax = 




3A2 (1 ) 

+ Sal) + [l40|/i3|^ + 15|?/t|^ — 167r(4ai + 15a2 + 200 : 3 )] > , 

A 


A 

167r2 

4|/^3p 

IGvr^ 


327r4 


\h3\\ 

A, 

A 2 


According to the SEHA2 convention [58] that we employ here. At is related to the “up” trilinear 
coupling {Tu )33 as At = {Tu) 33 /yt Note also that we choose A > 0, so that both A\ and nig 
are bound to assume exclusively negative values. 

At the leading order (one loop) gaugino masses at My are directly related to the gauge cou¬ 
pling constants and given by 

Mi = — fori = 1,2,3. (A.2) 

47r 


^Recall that intergenerational up trilinear couplings enter the scalar potential, F, as E D —HuQi{Tu)ij {XrIj . 
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As was shown in the numerical analysis of Sec. 3, some of the soft masses become negative 
for some values. As an example, one can derive from Eqs. (A.l) some naive bounds, where for 
simplicity all a’s are set equal to 1/20 at My, yt ~ 0.6, and tan /3 < 10 : 

^ or ^ < /i 3 < 1.3 . (A.3) 

However, the above values are just indicative. The full numerical bounds are shown in Sec. 3. 

B Sources of effective interactions below the messenger scale 

We will discuss here corrections to the dynamics of S due to various processes that might take 
place in a UV completion of our PQ NMMSM.^ We base our arguments on certain F-theoretic 
constructions. As explained, e.g., in Refs. [4, 5], chiral matter originates from two D7-branes 
intersecting along a 2-dimensional Riemann surface; matter interactions come from the intersection 
of three D7-branes at one point, which results in the intersection of three matter Riemann surfaces 
at this point. These intersections are generic, in the sense that they naturally emerge in the geometry 
of the string compactification. Thus, at the compactification scale the only natural matter couplings 
are cubic, and for this reason we set the mass term AiS' to zero in Sec. 2.2. The same argument 
also suggests that all natural couplings in the superpotential should be of the order of 1. The simple 
picture just described is generally amended in the presence of family effects [7], instantons, which 
we discuss below, and possibly other effects that we do not discuss here and are difficult to estimate 
without a detailed knowledge of the full F-theory model. 

We include in the analysis all important terms in the spurion X, whose vev gives a mass to 
the messengers and breaks SUSY. Below the messenger scale. My = {X), the spurion will be 
a nondynamical field, i.e., we set X —)■ {X) + 9'^Fx ■ Fx will be the primary source of SUSY 
breaking, which means that F-terms of other fields must be much smaller Fx in what 

follows). Note that {X) breaks also C/(1)pq. This means that at low energies there appear new 
effective f7(l)pQ breaking interaction terms. 

The contributions that might be relevant below My can have various sources. For the low- 
energy dynamics the most important will be terms with powers of the spurion superfield X, because 
X receives a large vev. As usual, operators of dimension 5 and higher will be suppressed by a large 
mass scale of the order of Mgut or even Mpp 

Giudice-Masiero term. Besides the superpotential term (2.4), the only cubic couplings of 5, X 
allowed by [7(1 )pq is in the Kahler potential: 

Yt 

5K = X^— - HuHd+ h.c., (B.l) 

Mgut 

where is a complex constant difficult to estimate [60]. 

The Kahler potential 6K can generate an additional y, term in the superpotential when the spu¬ 
rion X gets its vev, like in the Giudice-Masiero [41] mechanism: /tgm = A^F^/Mgut ■ However, 
one can either choose = 0, or perform a redefinition of the field, 5 —)• S — ^gm/A, which makes 

*For a different approach see, e.g., the analysis of Ref. [59]. 
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//gm disappear from the superpotential, and at the same time generates corrections to the tadpoles, 
and < 5 ^ 5 , and a Bfj, term, m 3 , of the form: 

5^s = -^rn'sAx - ^m|Ax, 

ml = -X^AxAx, (B.2) 

where m'g, m|, and Ax are given in Eqs. (2.9) and (A. 1). 

Note that in both cases the only origin of the /r-term remains the vev of the singlet, s = (5), 
so that iJefi = As. 

Install tons. In stringy models there are several types of instantons [61, 62]. An instanton action, 
Sx, depends on various moduli fields, which generally belong to two main categories, “twisted” 
and “untwisted”. Twisted instantons can transform by a shift under some abelian gauge symmetry, 
which in our case means that our U{1)pq effectively gives a charge to e~^^. When the moduli are 
set to some fixed value, becomes a numerical facfor, which we denofe here as Iq^n- In this 
case q is the C/( 1 )pq charge carried by the instanton and n indicates how many fields are associafed 
wifh if, as shown below. Some of fhe insfanfon confribufions are expecfed fo vanish, while fhe ones 
fhaf are presenf are expecfed fo be sfrongly suppressed, Iq^n ^ 1 - 

Insfanfon effecfs allow one fo generafe new ferms in fhe superpofenfial (corrections fo fhe 
Kahler pofenfial are negligible): 

5Wx = S + Mpi(/o ,2 xs + 1-2,2 s^) + Ii, 3 X^S + I-i, 3 XS‘^ + 1 - 3 , 3 ^^ + - (B.3) 

where Mpi is fhe reduced Planck mass and fhe ellipsis denofes higher powers of fhe fields, and we 
have here suppressed all ferms depending on X solely. 

As one can see, «: « 0, whereas ofher possible confribufions mighf be imporfanf due fo fhe 
large values of Mpi, {X), and Fx- The surviving ferms are 

5Wx A (MpV- 1,1 + Mpi(X)/o ,2 + h,3{X)^) S + (Mpi/_ 2,2 + I-1,3{X)) 52 

= Ml + r//o ,2 + ri^h,3) s + Mpi {1-2,2 + ti^- 1 , 3 ) , (B.4) 

5Lx = FxMpi (Jo ,2 + 2r?/i,3) S + I-i,3FxS^ , (B.5) 

where 77 = {X) /Mp\ ^ 1. 

From Eq. (B.4) follows fhaf if all Iq^n are approximafely of fhe same order of magnifude fhe 
mosf imporfanf confribufions are fhose characferized by fhe lowesf f7(l)pQ charges, q. However, 
if is nof uncommon for fhe dynamics of fhe UV sfringy fheory fo generafe selection rules fhaf can 
forbid certain ferms. If, for example, Iq<o,n = 0 due fo fhe UV dynamics, one is leff only wifh 

6Wx = {Xfh,3S, 

SLx = 2Fx{X)h,sS. (B.6) 

In our sfudy we always make fhis assumpfion, as if will be beneficial fo reducing fhe fine funing of 
fhe model. 
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Exchange of a heavy chiral 17(1 )pq neutral superfield. In general, string theory constructions 
imply the presence of some f7(l)pQ heavy chiral neutral superfields at the GUT scale that can be 
integrated out. 

The simplest superpotential coupling with one heavy chiral state 'h reads: W D Xq/XS'i! + 
Integrating out 'h gives 

which after SUSY breaking yields 

5W^ = —, (B. 8 ) 

and 

SLw,^ = ^Fx{X)S^ + h. c. = - Jrn'IS^ + h. c. (B.9) 

from the superpotential and 

SLk,^ = |5|2 (B.IO) 

from the Kahler potential (see Sec. 2.3). 

As was explained in Sec. 2.3, we expect ss Mqut and Aip « 1 if nonvanishing. 

Exchange of a 17(1 )pq heavy gauge hoson. Exchange of a heavy gauge boson V associated 
with the breaking of [/(1)pq can also contribute to the scalar Lagrangian: 

SKv = => SLv = forallj, (B.ll) 

where <I>j are the ?7(1)pq charged fields, qj are fheir charges, gpg is a generic coupling consfanf, 
and My ~ Mqut- Equafion (B.ll) implies fhaf all fhe sofl masses squared inj will receive an 
addifional correcfion of fhe order of gpQQxqjA.'^x'^ . 

Other contrihutions. There are more terms preserving 17(1)pq which mighf renormalize fhe 
Kahler pofenfial as, for example, 6K = ^X^XXS + c. c. Wifh M « Mqut they are negligible, 
so fhaf we will nol fry fo pinpoinf fheir sources here. 


Leading supergravity corrections. Supergravify corrections [63] for W = VTo + FxX + ^fS + 
2 /r'5'^ (with ^F <C Fx) and for the canonical Kahler potential subject to the condition of vanishing 
cosmological constant read: 


()U(5) = 



\/3Afpi 



+ qliJkL 15-12 
^ 3M2 


+ .. 


(B.12) 


They are also negligible. 
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In summary, when expressed in terms of A = Fx/My and x = My the full correc¬ 

tions to the superpotential and soft SUSY breaking terms are the following: 


5mf 

f 


m'l 

ml 


ip 


is 


gpQqxqi\J^?x^ 

-lA^PlAl^x^ + g%Q qxqs\M'^x^ 

—A|,Mgut 

2A^AAfGUT = —2g!A , 

-X^AxAx , 

-2/uAMa„T^" - ^A2a:=>MouT + ^A“x=> 


-^rnlAx, 


(B.13) 


where in the last equation of (B.13) m| is the term calculated in Appendix A. 

Since typical values of A imply that Ax <C MguTj in our study we assume A^ = 0 without 
loss of generality. Moreover, A^x^ <C Afl^gy ^o that we can neglect the gpg corrections to the soft 
masses. Note that when A^ = 0 and gpq = 0, Eqs. (B.13) reduce to the form given in Eqs. (2.9). 


C Estimates of constants 


In this appendix we estimate the order of magnitude of different coefficients appearing in the Ea- 
grangian. Eor simplicity in what follows all coefficients are assumed to be real. We consider con¬ 
straints that come from the tree-level vacuum equations of motion (e.o.m.), and those that follow 
from the stability of the physical vacuum. 

Eet us recall that k = 0, = 0, and we are interested in small tan /3 values. Moreover, we 

assume that all the GMSB soft masses given in Appendix A are of the order of Msusy ^ A. We 
also assume that no accidental cancellation takes place between different terms. We only consider 
terms in leading powers of v, and do not explicitly differentiate between constants of the order of a 
few. 

The vacuum e.o.m. are given by: 


2(As Ax + ml) 


sin 2/3 - ~ ^ = 0 . 

2X^s^ + + mir +m% 

Fid Flu 

[X^-g"^) -ptan2/3 -f (As Aa- f m|) = 0, 


tan 2/3 


s (m| -f m'l -f /r'2 -f A^ v^) -h - Xv^ (Aa + g!) = 0 . 


(C.l) 

(C.2) 

(C.3) 


where = {gf -|- < 72 )/^ ^3 = + X{g's + ip). Short inspection of the e.o.m. reveals 

that for tan /3 equals a few, s is of the order of the soft SUSY masses, s ~ Msusy- Thus, since 
sin 2/3 ~ 1 one also finds m| ~ g's + Cf ^ Msusy- 

On fhe ofher hand, fhe requiremenf of free-level sfabilify in fhe scalar secfor yields ^ 


[-2»A + {A,+i,') sin 2/3] N - (A±M1. 


(C-4) 


’We analyze here the diagonal minors of the scalar mass matrix, which can be explicitly found, e.g., in Ref. [40]. 
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which implies 


“ {^s + iplj!)/s > max(/i'^, M|usy) ■ (C.5) 

By using Eq. (C.3) one can derive 


(m| + m's + > max(/i'^, M|usy) > 


(C. 6 ) 


which is trivially satisfied when m'g > 0 and rrig < 0, as is the case at hand (see Appendix A). On 
the other hand, requiring tree-level stability of the pseudo-scalar masses, whose expressions can 
also be found in [40], yields an additional constraint: 


— 2m'g — 




> 0 . 


(C.7) 

Thus Eq. (C.3) also implies (m| — m'J -|- /r'^) > 0, which in turns implies m'g < . Assuming 

a relation m'J = —c/r'A, with c > 0, just like in Eq. (2.9), one gets 




(C. 8 ) 


These arguments show that for small tan (3, which is the case of interest here, one can get 
reasonable physics only when the mass scale associated with ^ 5 , m'J, and /r' is much larger 
than Msusy- 

By reinserting Eq. (C. 8 ) back into Eqs. (C.1)-(C.3) one can finally infer that s + Cf/ fr' ~ 0 
and s + Cs/fT's ~ 0 should simultaneously hold, up to terms that are much smaller than Msusy- 
Thus, one can derive the relation (2.11), 


Cs _ m'l 
ip 


(C.9) 


which must hold at the EW scale. 

To summarize, we present all the coefficients for which the hierarchy Msusy < A ^ 1 /^' 
holds in terms of Msusy, A, and |//'|: 


~ -A^suSYf^^ is ~ MsusyAfi', -A/i' (C.IO) 


Our numerical analysis confirms these predictions. Note that ^' < 0 implies ip > 0. The sign 
of ip (as well as that of A) can be always be chosen positive by rescaling the S and, e.g., Hu 
superfields. 

As we claimed in Sec. 2.3, the relation (C.9) is set at the My scale when appropriate instanton 
corrections are chosen. More importantly for the fine funing calculation, however, Eq. (C.9) is 
approximately RG invariant, so that it continues to hold at different scales. We prove this by 
writing down the one-loop leading terms in the RG running: 

51og^i;’ = A^, (51og/i' = 2A^ 

Slogis ~ A2(1 -h 2AxiF/is), (51ogm'| = 2X‘^{1 -f 2Axn'Im '^), 

where we have used the estimates in (C.IO) and the simplified notation <5 = IGvr^^. The terms in 
parentheses are approximately 1 so that 

~ <^log PS 0, (C.ll) 

which leads to Eq. (C.9). 
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D Tree level Higgs mass in the PQ NMSSM 


Tree-level masses of the three CP-even Higgs bosons in the PQ NMSSM can be calculated by di¬ 
agonalizing the mass matrix explicitly given, e.g., in Ref. [40]. We give here approximate formulas 
for the lightest eigenvalue. 

We assume the following hierarchy of the mass matrix eigenvalues: <C rrih^ <C ruh^. 

One can thus show that 


m 


Til,tree 


2 / 2 2 no , ^2 ■ 2 no^ , (2sA - sin 2^ -//'sin 2/3)" 

7 ( g cos 2/3 -h A sm 2/3)4- —-■ ■ - --—— 

^S + 


(D.l) 


Since g! is much larger than the GMSB-generated soft terms given in Eqs. (A.l), one can 
further reduce Eq. (D.l) to 


^ hi,tree 


Mf cos^ 2/3 4- A^u^ sin^ 2/3 





12 


{is + iF^J^')/s 


(D.2) 


Note that in the model presented in this paper (^5 + iF^J‘')/s < 0. 

One can schematically separate Eq. (D.2) into the MSSM and A-dependent parts: 

/Vf| cos^ 2/3 4- . One should recall from Appendix C that —g's, so that we get 


(5m? , = A^u^ sin^ 2/3 f 1-^— | Rs sin^ 2/3 • , (D.3) 

( 1 + ^/ jp,.' 

which was used in the discussion of Sec. 3. As was mentioned over there, when ^ one 
gets (5m^^ « 0, even for large values of A > 0.7, which is the case presented in BP2 of Table 1. 

Recall, finally, from Appendices B and C that —^s = 2A^i?, and 2A < —g ', which implies 
that in our model e = is/iFfJ-' ^ 1- As we discussed in Sec. 3, we found no points for which 
m|i > /\T| , so that large radiative corrections to the Higgs mass are always necessary to 
obtain nihi ~ GeV. 
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